1 In the claims: 

2 1 . A method for gamut mapping of an input image using a space varying 

3 algorithm, comprising: 

4 receiving the input image; 

5 converting color representations of an image pixel set to produce a corresponding 

6 electrical values set; 

7 applying the space varying algorithm to the electrical values set to produce a 

8 color-mapped value set; and 

9 reconverting the color-mapped value set to an output image. 

1 0 2. The method of claim 1 , wherein the space varying algorithm minimizes a 

1 1 variational problem represented by: 

12 E(u) = j"(-D 2 + a\VD\ 2 )dQ. , subject toueS, wherein Q is a support of the input 

a 

13 image, a is a non-negative real number, D=g*(u-u 0 ), g is a normalized Gaussian kernel 

14 with zero mean and a small variance a, Uq is the input image, and u is the output image. 

15 3 . The method of claim 2, further comprising: 

1 6 solving the variational problem at a high value of a; 

17 solving the variational problem at a low value of a; and 

18 averaging the solutions. 

19 4. The method of claim 3 , wherein the step of averaging the solutions comprises 

20 using a spatially adaptive weighting scheme, comprising: 

u fmal [k, j] = w[k, j]u small [k, i](l - w[k, y [*> A 

21 wherein the weight w[k, j], comprises: 

w\k, j]= — ^ a ' and 

l + /l\Vg*u 0 \ 2 

22 wherein p is a non-negative real number. 

23 5. The method of claim 2, wherein the variational problem is solved according to: 

24 — = ag * AD - g * D, subject to u e <9 . 
dt 

25 6. The method of claim 2, wherein the space varying algorithm is solved according 

26 to: 
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1 = u "j + t[ccL"j - D" ), subject to u" } e 3, wherein 

r = dt, 

^ D" = g*g*(a"-« 0 ) 

i" = £> 2 * (w" - u 0 ) and 

A = 

3 7. The method of claim 1 , wherein the space varying algorithm minimizes a 

4 variational problem represented by: 

5 E(u) = j(pi(D)+ ap 2 (\VD\))dQ. , subject to u e &, wherein p x and p 2 are scalar 

a 

6 functions. 

7 8. The method of claim 2, further comprising: 

8 decimating the input image to create one or more resolution layers, wherein the 

9 one or more resolution layers comprises an image pyramid; and 

I o solving the variational problem for each of the one or more resolution layers. 

II 9. The method of claim 1 , wherein the method is executed in a camera. 

12 10. The method of claim 1 , wherein the method is executed in a printer. 

13 11. A method for color gamut mapping, comprising: 

14 converting first colorimetric values of an input image to second colorimetric 

1 5 values of an output device, wherein output values are constrained within a gamut of the 

16 output device; and 

17 using a space varying algorithm that solves an image difference problem. 

18 12. A computer-readable memory for color gamut mapping, comprising an instruction 

19 set for executing color gamut mapping steps, the steps, comprising: 

20 converting first colorimetric values of an original image to second colorimetric 

21 values, wherein output values are constrained within a gamut of the output device; using a 

22 space varying algorithm that solves an image difference problem; and 

23 optimizing a solution to the image difference problem. 

24 13. The computer-readable memory of claim 12, wherein the image difference 

25 problem is represented by: 

26 E(u)= \{D 2 +a\VD\ 2 )dCl 

n 



HP 10005732 



-13- 



1 subject to u e &, wherein Q. is a support of an input image, a is a non-negative real 

2 number, D=g*(u-Uo), g is a normalized Gaussian kernel with zero mean and small 

3 variance a, Uo is the input image, and u is an output image. 

4 14. The computer-readable memory of claim 12, wherein the instruction set further 

5 comprises steps for: 

6 solving the image difference problem at a high value of a; 

7 solving the image difference problem at a low value of a; and 

8 averaging the solutions. 

9 15. The computer-readable memory of claim 14, wherein averaging the solutions 

10 comprises using a spatially adaptive weighting scheme, comprising: 

«jw [*, 7'] = w[k, j]u small [k, j](l - w[k, j])u klgh [k, j\ and 

1 1 wherein the weight w[k, j], comprises: 

12 w[£,y']= — ^ ^,and 

l + /]\Vg*u 0 \ 2 

13 wherein p is a non-negative real number. 

14 16. The computer-readable memory of claim 1 2, wherein the image difference 

1 5 problem is represented by: 

16 E(u)= ^(0)+ ap 2 (\VD\))dQ., wherein p x and p 2 are scalar functions. 

n 

17 17. The computer-readable memory of claim 12, wherein the instruction set further 

18 comprises steps for: 

19 decimating the input image to create one or more resolution layers, wherein the 

20 one or more resolution layers comprise an image pyramid; and 

2 1 solving the image difference problem for each of the one or more resolution 

22 layers. 

23 18. The computer-readable memory of claim 1 7, wherein the instruction set further 

24 comprises steps for: 

25 (a) initializing a first resolution layer; 

26 (b) calculating a gradient G for the resolution layer, the gradient G comprising: 



HP 10005732 



-14- 



G = A(u -u 0 )+a K (u-u o ), wherein Ax is a convolution of each color 



2 plane of x with K L . 



0 1 0 

1 -2 1 
0 1 0 



and ak=a 0 *2 2(k " 1) ; 



3 



4 [i 0 is a constant; 



(c) calculating a normalized steepest descent value L,- Lj.i-|a 0 VNSD*G, wherein 



5 (d) projecting the value onto constraints Proja (Lj), wherein Proj 9 (x)is a projection 

6 of x into a gamut 3; and 

7 (e) for a subsequent resolution layer, repeating steps (b) - (d). 

8 19. A method for image enhancement using gamut mapping, comprising: 

9 receiving a input image; 

1 0 from the input image, constructing an image pyramid having a plurality of 

1 1 resolution layers; 

12 processing each resolution layer, wherein the processing includes completing a 

13 gradient iteration, by: 

14 calculating a gradient G; 

15 completing a gradient descent iteration; and 

1 6 projecting the completed gradient descent iteration onto contstraints; and 

17 computing an output image using the processed resolution layers. 

1 8 20. The method of claim 19, wherein the gradient G, comprises: 

19 G = A(u-u 0 ) + a k (u-u 0 ), 

20 wherein u is the output image, u 0 is the input image, and a is a non-negative real 

21 number. 

22 21. The method of claim 1 9, wherein completing the gradient descent iteration 

23 comprises calculating: 

24 = SG /i(C7* AG) + a k ZG>) > ^ 

25 Lj =L j _ 1 -<u 0 -Mnsd- g > 

26 wherein ju NSD is a normalized steepest descent parameter, ju 0 is a constant, k is a number 

27 of resolution layers in the image pyramid, and j is a specific resolution layer. 

28 22. The method of claim 1 9, wherein proj ecting the completed gradient descent 

29 iteration onto the constraints is given by: 
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1 L J =Proj a (L J ), 

2 wherein Proja(x) is a projection of x into a gamut 

3 23. The method of claim 19, wherein constructing the image pyramid, comprises: 

4 smoothing the input image with a Gaussian kernel; 

5 decimating the input image; and 

6 setting initial conductive L D = max{Sp}, wherein Sp is an image with the coarsest 

7 resolution layer for the image pyramid. 

8 24. The method of claim 23, wherein the Gaussian kernel, comprises: 
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10 25. The method of claim 19, wherein processing each resolution layer further 

1 1 comprises applying a space varying algorithm to minimize a variational problem 

12 represented by: 

13 E(u)= j(l) 2 +a\vD\ 2 )dQ, subject to wherein Q is a support of the 

n 

14 image, and D=g*(u-Uo), wherein g is a normalized Gaussian kernel with zero mean and 

1 5 small variance a, u<, is the input image, u is the output image, and wherein a is a non- 

1 6 negative real number. 

1 7 26. The method of claim 1 9, wherein processing each resolution layer comprises 

1 8 applying a space varying algorithm to minimize a variational problem represented by: 

1 9 E(u) = (D) + ap 2 (\VD\))dQ. , subject to u e S , wherein p l and p 2 

a 

20 are scalar functions. 

21 27. The method of claim 26, wherein p x and p 2 are chosen from the group 

22 comprising p(x) = |x| and pix)- yjl + x 2 . 
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